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Abstract 

A connection between fractional calculus and statistical distribution theory 
has been established by the authors recently. Some extensions of the results 
to matrix- variate functions were also considered. In the present article, more 
results on matrix-variate statistical densities and their connections to frac- 
tional calculus will be established. When considering solutions of fractional 
differential equations, Mittag-Lefner functions and Fox H- function appear 
naturally. Some results connected with generalized Mittag-Lemer density 
and their asymptotic behavior will be considered. Reference is made to 
applications in physics, particularly superstatistics and nonextensive statis- 
tical mechanics. 
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1. Introduction 

For our discussions to follow we need the definitions of Mittag-Lefner 
function, Wright function, Fox H-function and generalized Mittag-Lemer 
density. Hence we will introduce these functions first. 

The Mittag-Lefner function is used in many areas of physical sciences 
[see Haubold, Mathai, and Saxena (2009), Saxena, Mathai, and Haubold 
(2010), Haubold, Mathai, and Saxena (2011)]. Haubold and Mathai [2000, 
see also Haubold and Kumar (2010)] have shown that when we move from 
total differential or integral equation to fractional cases the exponential 
type solutions go in terms of Mittag-Lemer functions. For a modern intro- 
duction to fractional calculus see Mainardi (2010) and Tenreiro, Machado, 
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Kiryakova, and Mainardi (2011). Solutions of certain fractional reaction- 
diffusion equations belong to Mittag-Leffler functions [Gorenflo and Mainardi 
(1996), Kiryakova (2000), Kilbas (2005), Mainardi and Pagnini (2008), 
Mathai, Saxena, and Haubold (2010)]. Let us start with the definition of 
Mittag-Leffler function as proposed by Mittag-Leffler and its various types 
of generalizations. 

1.1. Mittag-Leffler Function 

The basic Mittag-Leffler function is an extension of the exponential func- 
tion. The exponential series is given by the following: 



If we make a change in T(l + k) and write it as T(l + ak) for 3?(a) > 0, 
where means the real part of (•) then we have the original Mittag-Leffler 
function, namely, 



Thus for a = 1 we have E\{z) = e z . Three forms of Mittag-Leffler functions 
are frequently used in the literature and hence we list these three forms first. 



E ^ {x) = £ k\r(p + ak) ' (7)fc = 7(7 + 1} - (7 + k " 1} ' (7)o = 1,7 * °' a > °' p > 0; 
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When the parameters are in the complex domain the conditions become 
3ft(a) > 0, > 0. It is easy to note that the Mittag-Leffler functions 

of (l)-(4) are special cases of the Wright function [Wright (1940)] which is 
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defined as 



P^ q {z) = P ipq(ai, .., a p ; h, b q ; z) 

n; = i r(o 3 - + k) z k 



where the ctj's and fy's are complex parameters. This Wright function is a 
special case of the most generalized special function, namely, the H-function 
and the H-function is defined by the following Mellin-Barnes integral: 

r I \ I M 1 fc+ioo 

H m ' n (z)=H m ' n z \y ai > ai >>->( a P' a p) = _ / 6(s)z~ s ds (6) 



where 



inr = i nbj + Pjsmuu r (* - ^ - M> 

{ILUh-i r(i - 6, - MHUU+i T ^ + 



where ctj's and 6j's are complex quantities, ay's and /3j's are positive real 
numbers, i = y 7 — 1. Various types of possible contours, convergence con- 
ditions and applications may be seen from Mathai, Saxena and Haubold 
(2010). The elementary special functions such as binomial function, Bessel 
functions of various types, sine and cosine functions, hyper geometric func- 
tions of different types etc, which are applicable in physical sciences, are 
special cases of (6). 

1.2. A Generalized Mittag-Leffler Statistical Density 

A statistical density in terms of Mittag-Leffler function was originally 
defined by Pillai (1990) in terms of the following distribution function or 
cumulative density: 

G y (y) = 1 - E a (-y a ) = p ^-L- J fc) , < a < 1, y > (8) 

and G y (y) = for y < 0. Since it is differentiable one can obtain the density 
function by differentiating on both sides with respect to y and the density 
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function f(y) is the following: 

f(y) = l y Gy{y) 



d_ 

dy 



^ r ( 1 + ^) 



y, (-l) k+1 aky ak - 1 



by replacing k by k + 1 

/(y) = y^E^i-y"), < a < 1, y > (10) 

where E a ^(x) is the generalized Mittag-Leffler function. 

In order to study the properties of the Mittag-Leffler function, the cor- 
responding fractional differential or integral equations and their generaliza- 
tions, usually the techniques from Laplace and Fourier transforms are used. 
Hence we will consider some Laplace transforms and then we will establish 
some connections to Levy distributions. The Laplace transform of (10) is 

/>00 /'OO 

L f (t)= / e- tx f(x)dx= / e- tx x a - l E a , a {-x a )dx (11) 
Jo Jo 

= (i + t a y\\e\ < 1. (12) 

Note that (12) is a special case of the general class of Laplace transforms 
discussed in Section 2.3.7 of Mathai and Haubold (2008). From (12) one 
can also note that there is a structural representation in terms of positive 
Levy distribution. 

Definition 1.1. A positive Levy random variable u > 0, with parameter 
a is such that the Laplace transform of the density of u > is given by 
e~ r . That is, 

E[e~ tu } = e~ ta (13) 
where E(-) denotes statistical expectation or the expected value of (•). 

When a = 1 the random variable is degenerate with the density/ prob- 
ability function 

J 1, for ar = 1 

fi( x ) = { ( 14 ) 
0, elsewhere. 
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Consider a gamma random variable with the scale parameter 5 and shape 
parameter (3 or with the density function 



x l3-l e -S 



/l(x) = j wm ' for - x < > M > (15) 

[0, elsewhere; 
and with the Laplace transform 

L fl (t) = (l + 6t)-e. (16) 

A more general form of the Mittag-Lefher density is given by the fol- 
lowing [Mathai (2010)]: 



^ fc!r(a?? + afc) a^+ fc 1 J 

™ar?-l / „a\ 

— E ^{--) < 18 > 

and with the Laplace transform 

L g (t) = [1 + a^H, \at a \ < 1 (19) 

where a > 0, < a < 1, 77 > 0. 

With these preliminaries we will discuss a few problems connected with 
Mittag-Lemer density and we will examine the matrix-variate analogues. 

2. Generalized Mittag-Leffler Density 

First, we will examine a minor generalization of a known result, which 
is given recently by Mathai (2010). The proof is given by others by using 
different methods but we will obtain it by using some properties of condi- 
tional expectations. To this end, we will list a basic result on conditional 
expectations in the form of the following lemma. 

Lemma 2.1. For two random variables u and v having a joint distribution, 

E{u) = E[E(u\v)] (20) 

whenever all the expected values exist, where the inside expectation is taken 
in the conditional space of u given v and the outside expectation is taken in 
the marginal space of v. 
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A real scalar random variable x is a mathematical variable where a 
probability statement of the type Pr{x < a}, that is the probability of 
x falling in the interval — oo < x < a, is defined for all real values of a. 
A density function f(x), associated with the random variable x, is a non- 
negative integrable function with the total integral unity. The expected 
value of a function <j>{x) of x, denoted by E[<f>(x)], is the integral 

/oo 
4>{x)f(x)dx (21) 
-oo 

where f(x) is the density function of x when x has nonzero probabilities on 
a continuum of points. Parallel definitions go for discrete random variables, 
joint densities of many random variables, vector or matrix random variables 
etc. 

The lemma follows from the definition itself. By using this lemma and 
equations (13), (15) and (17) one can obtain a structural representation for 
a Mittag-Lefner random variable. This will be stated as a theorem. 

Theorem 2.1. Let y > be a Levy random variable with Laplace transform 
as in (13), x > be a gamma random variable with the density as in (15) 
and let x and y be independently distributed. Then u = yx& is distributed 
as a Mittag-Leffler random variable with Laplace transform 

L u (t) = [l + 8t a ]- p . (22) 



Proof: For proving this result we will make use of the conditional argument 
as given in Lemma 2.1. Let the density of u be denoted by g(u). Then the 
conditional Laplace transform of g, given x, is given by 

E[ e -^y\x] =e~ tax . (23) 

But the right side of (23) is in the form of a Laplace transform of the density 
of x with parameter t a . Hence the expected value of the right side, with 
respect to x, is available by replacing the Laplace parameter t in (16) by t a , 

L g (t) = (l + 5t a )-P,5>0,l3>0 (24) 

which establishes the result. From (20) one property is obvious. Suppose 
that we consider an arbitrary random variable y with the Laplace transform 
of the form 

L g (t) = e-W*)] (25) 
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whenever the expected value exists, where 4>(t) be such that 

<f>(txa) = x<f>{t), Jim <f>(t) = 0. (26) 
Then from (23) we have 

E[e-<- tx °')v\ x ] =e -*W)] m (27) 

Now, let x be an arbitrary positive random variable having Laplace trans- 
form, denoted by L x (t) where L x (t) = tp(t). Then from (25) and (20) we 
have 

L g (t) = i/>W)]- (28) 

General properties of Mittag-Leffler functions, their generalizations, their 
applications to fractional calculus, reaction-diffusion type problems, may 
be seen from papers, for example, Gorenflo and Mainardi (1996), Gorenflo, 
Kilbas, and Rogosin (1998), Gorenflo, Loutchko, and Luchko (2002), Kilbas 
(2005), and Kilbas and Saigo (1996). 

We can create a generalized Mittag-Leffler density as follows: Consider 
the generalized Mittag-Leffler function 



01 ( x ) - JL V (~ 1 ) fcr ( ? ? + fc ) ar>-l+afc (2g) 



x a "- 1 ^ J)a (-x a ),a>0,»/>0. (30) 



Then the Laplace transform of gi(x) is the following: 



k=0 

This is a special case of a general class of Laplace transforms considered in 
Mathai, Saxena, and Haubold (2006), Saxena, Mathai, and Haubold (2010), 
and Haubold, Mathai, and Saxena (2011). 



2.1. Mellin-Barnes Representations 
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It is shown in Mathai (2010) that for handling problems connected with 
Mittag-Lemer densities it is convenient to use the Mellin-Barnes represen- 
tation of a Mittag-Leffler function and then proceed from there. gi(x) of 
(30) can be written as a Mellin-Barnes integral and then as an H-function. 
Then 



9i(x) 



1 1 

T(r?) 2jH 

ari-l 

n \,1 



C+JOO 



T(s)T(r]-s] 



T(rf) " il2 r l(0,l),(l-a»7,a 

= _2 l r 1+100 nv-h+znh-z) 

aT(rj) 2tti / cl _ ioo 
by taking arj — 1 — as = —s±, 



x ai i- l (x a )~ s ds, < c < sR(r?) (33) 

(34) 

x~ s ds (35) 



T(l - s) 



9i(x) 



i.i 



la- 1 , 1 ) 

1(^-1 ,1),(0,1) 



(36) 



One can also look upon the Mellin transform of the density of a positive 
random variable as the (s — l)-th moment and therefore, from (2.15) one 
can write the Mellin transform as an expected value. That is, 



M gi {s) = E(x s - 1 ) in 9l 

_ v ' a. a' y a a' 

T(ri) ' aT(l -s) ' 

T( V ) r(2 - s) 

for 1 — a < 3?(s) < 1, < a < 1, 77 > 0. Observe that the right side in (39) 
goes to 1 when s — > 1 thereby establishing that the non-negative function 
gi(x) is a density function. This is the most convenient way of showing 
that gi(x) is a density because term by term integration from to 00 is not 
possible in a Mittag-Lemer function, which is in the form of a series. This 
gi (x) is called the generalized Mittag-Lemer density. 

From (36) we may also observe the following properties. If u is the 
generalized Mittag-Lemer variable then from (35) and (36) we have 



(37) 
(38) 

(39) 
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E(u- 



s-1 



) = 



i rfo-i + a)r(i + i-a) 
r(r?) r(2 - S ) 



(40) 



= J B(y s - 1 ) J E(x-) 



s-1 



(41) 



where 




-) 



(42) 



and 



Eix-y- 1 = 



(43) 



which suggests that the Mittag-Leffler random variable has the representa- 
tion 



where y is a positive Levy random variable with the (s — l)-th moment in 
(42), x is a gamma random variable with the density function in (15), u 
is a Mittag-Leffler variable and it is assumed that x and y are statistically 
independently distributed. From (42) it may be further noted that y does 
not have a density but for large values of y one can have an asymptotic form 
of the density of y. 

3. The Pathway Model 

In model building situations in physical sciences, what is usually done 
is to select a model from a parametric family of distributions if it is a 
deterministic (non-random) or random (non-deterministic) situation, such 
as a gamma family of distributions. A gamma family of distributions usually 
has two parameters (a, /3). Thus these two parameters are suitably chosen 
and then a model is selected for the data at hand. Often it is found that 
the selected model is not a good fit because the data requires a model with 
a thicker or thinner tail than the one given by the parametric family, or it 
may be that the appropriate model is in between two parametric families 
of distributions. In order to come up with an appropriate model in such 
situations a pathway model was proposed in Mathai [2005, see also Mathai 
and Haubold (2010), Mathai, Haubold, and Tsallis (2010)]. The original 
pathway model of Mathai (2005) is for the rectangular matrix cases. The 
scalar version of the pathway model is the following: 



u = y(x<* 



) 



(44) 
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f 1 (x) = c 1 \x\"<[l-a(l-a)\x\ s }^ (45) 

for a > 0, r] > 0, 5 > 0, 7 > — 1 and ci is a constant. A statistical density 
can be created out of fi(x) under the additional conditions l-a(l-a)\x\ s > 
and then c\ will act as a normalizing constant. Hence the model is directly 
applicable to deterministic situations as well as random or non-deterministic 
situations. For a < 1 the model in (45) will stay in the generalized type-1 
beta family of functions. For a > 1, writing 1 — a = —(a — 1) the model 
switches into the model 

h{x) = c 2 M 7 [l + a(a - l)\x\ 5 }'^ (46) 

for a > 1, —00 < x < 00, a > 0, rj > 0, S > 0. Note that when a switches 
from a < 1 to a > 1 the model goes from a generalized type-1 beta family 
of functions to a generalized type-2 beta family of functions. Again, taking 
C2 as the normalizing constant one can create a statistical density out of 
f2(x). When a goes to 1 from the right or from the left, f\{x) and f2(x) 
will go into a generalized gamma family of functions, namely, 

lim h{x) = lim f 2 (x) = f 3 (x) = cslxPe-^l* (47) 

where a > 0, r\ > 0, 5 > 0. It is shown in Mathai and Haubold (2007) 
that almost all statistical densities in current use in statistics and physics 
are special cases of the pathway model of (45) to (47) or as a one-to-one 
function of the variable therein. The matrix-variate pathway model is of 
the following form: 

f(X) = c\A^{X - M)B(X - M)'A^[< 

x I J - a(l - a)A2 (X - M)B{X - M)'A* | ^ (48) 

where a prime denotes a transpose, M, A, B are constant matrices, M is 
m x n, m < n, A = A' >0 is m x m and real positive definite, B = B' > 
isnxn and real positive definite, X is mxn matrix of distinct real variables 
and of rank m, and |(-)| denotes the determinant of (•). Note that for m = 1 
one has a quadratic form in pathway vectors and then (48) can produce an 
extended density of quadratic forms, see Mathai (2007). 

One can establish a path of going from a Mittag-Lefner variable to a 
positive Levy variable. Replace a by a(q — 1) and 77 by rj/(q — 1) in the 
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Laplace transform (22) of the generalized Mittag-Leffler density in (8) and 
(18). That is, 

L g (t) = [l + a(q-l)t a ]-^,q>l. (49) 

If q — > \ + then 

L g (t) ->• e- a "* Q = L u (t) (50) 

which is the Laplace transform of a constant multiple of a positive Levy 
variable with parameter a, given in (13). Thus q here creates a pathway of 
going from the general Mittag-Leffler density g to a positive Levy density 
with parameter a, the multiplying constant being (ar])a. For a discussion 
of a general rectangular matrix- variate pathway model see Mathai (2005). 
The result in (48) can be put in a more general setting also. These ideas are 
generalized to the matrix- variate cases recently, see Mathai (2010). Some 
multivariable, not matrix variable, analogues may be seen from Lim and 
Teo (2009) and Pakes (1998). 

Another area of applications of Mittag-Leffler functions, Mittag-Leffler 
density and pathway models is in the area of quadratic and bilinear forms 
in Gaussian variables. The theory and applications may be seen from the 
books of Mathai and Provost (1992) and Mathai, Provost and Hayakawa 
(1995). With the help of the pathway models of (45) to (47) one can extend 
the theory of quadratic and bilinear forms to wider classes, rather than 
confining to Gaussian variables. Some works in matrix variate cases in this 
direction may be seen from Mathai (2007). 

Particular cases of the models in (45) to (47) are Tsallis statistics and 
superstatistics of non-extensive statistical mechanics. In (45) and (46) put 
7 = 0, 5 = 1, a = 1, r] = 1 and consider only positive variable x > 
then we get the statistics of Tsallis (1988, 2004, 2009). Superstatistics is 
not available from the form in (45) for a < 1. Consider (46) for x > 0. Put 
a = 1, 5 = 1, 7] = 1 then we get the superstatistics of Beck and Cohen 
(2003) and Beck (2006, 2008, 2009, 2010). 

The models in (45) to (47) can also be derived by optimizing a gener- 
alized measure of entropy under some moment-type conditions, see Mathai 
and Haubold (2007). 

4. Multivariable Analogues 

A slight extension of the Levy variable is the Linnik variable. A Lin- 
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nik random variable is denned as that real scalar random variable whose 
characteristic function is given by 

= TTjtp 0<a ^ 2 ' -°°<t<™- (51) 

For a = 2, (51) corresponds to the characteristic function of a Laplace 
random variable and hence Pillai and his coworkers name the distribution 
corresponding to (51) as the a-Laplace distribution. Infinite divisibility, 
characterizations, other properties and related materials may be seen from 
the review paper Jayakumar and Suresh (2003) and the references therein, 
Pakes (1998) and Mainardi and Pagnini (2008). Multivariate generaliza- 
tion of Mittag-Leffier and Linnik distributions may be seen from Lim and 
Teo (2009). Since the steps for deriving results on Linnik distribution are 
parallel to those of the Mittag-Leffier variable, further discussion of Lin- 
nik distribution is omitted. The a-Laplace distribution plays vital roles in 
non-Gaussian stochastic processes and time series. 

A multivariate Linnik distribution can be defined in terms of a multi- 
variate Levy vector. Let T' = (ti, t p ), X' = (xi, x p ), prime denoting 
the transpose. A vector variable having positive Levy distribution is given 
by the characteristic function 

E[e iT ' x }=e^ T ' i:T ^, 0<a<2, (52) 

where S = S' > is a real positive definite p x p matrix. We consider a 
representation, corresponding to the one in (44). 

U = y«X (53) 

where the pxl vector X, having a multivariable Levy distribution with pa- 
rameter a, and y a real scalar gamma random variable with the parameters 
5 and /3, are independently distributed. Then the characteristic function of 
the random vector variable U is given by the following: 



E[e iy " T ' x ] = E[E[e iylaT ' x ]\ y ] (54) 

= E[e-y^ T ^} = [1 + <J(T'Sr)f (55) 

Then the distribution of U, with the characteristic function in (55) is called 
a vector-variable Linnik distribution. Some properties of this distribution 
are given in Lim and Teo (2009). 
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One can also establish some central limit property but it will be of 
the nature of generalized Mittag-Leffler variable going to a positive Levy 
variable. Consider the generalized Mittag-Leffler density with the Laplace 
transform 

L x (t) = [1 + 5t a ]-P, 5>0, p > 0, < a < 1. (56) 

Let Xj,j — 1, ...,n be independently and identically distributed as in (56) 
and let 

(xi + ... + x n ) 

y = t 

Then 

L y {t) = [l+°±-]-^. (58) 
Hence when n — > oo 

lim LJt) = e~ 5l3ta (59) 
n— k» 

which is the Laplace transform of a constant multiple of a positive Levy 
variable with parameter a, the constant being (<5/3)«. Thus the central 
limiting property is that a certain normalized sample mean from a Mittag- 
Leffler population goes to a Levy variable. One can also introduce the 
pathway model here. If we replace <5 by S(q — 1) and j3 by Q > 1 then 
as q — > 1 , 

lim[l + 5(q - l)t a Y^ = e-^* Q (60) 

q-tl 

where q is a pathway parameter describing the path of going from a general 
Mittag-Leffler variable to a constant multiple of a Levy variable. 

We can extend the above ideas to matrix-variate cases also. Let X = 
(x rs ) be m x n where all distinct, X be of full rank, and having a 

joint density f(X), where f(X) is a real- valued scalar function of X. Then 
the characteristic function of f(X), denoted by (f>x{T), is given by 

0x(T) = E[e^ XT \i = y^l and T = (t rs ) (61) 

is an n x m matrix of distinct parameters t rs 's and let T be of full rank. As 
an example, we can look at the real matrix-variate Gaussian distribution, 
given by the density 



(57) 
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g(X) = \3l^ e -*(AXBx>) (62) 

7T 2 

where X is m x n, A = A' > 0, B = B' > are m x m and n x n positive 
definite constant matrices, a prime denoting a transpose. Then the Fourier 
transform of g(X) is given by the following: 

<f> g (T) = E[e it < XT \i = v^T (63) 

= e -|tr(B-^TA- 1 T'B-^) ^ 

where T is the parameter matrix, and tr(-) denotes the trace of (•). Moti- 
vated by (61) and (64), Mathai (2010) defined matrix- variate Linnik density 
and gamma-Linnik density. We will list some of them here. 

Definition 4.1. Matrix- variate Linnik density. Let X and T be as 

defined in (61). Then X will be called a real rectangular matrix- variable 
Linnik random variable if its characteristic function is given by 

<f> x (T) = e-IMTEiT^)]^ < a < 2 (65) 

where Si > is m x m and £2 > is n x n positive definite constant 
matrices. 

For the real rectangular matrix- variate Gaussian density of (62), Si = 
\A~ X and S2 = \B~ X and a = 2. The following result is given recently in 
Mathai (2010). 

Theorem 4.1. Let y be a real scalar random variable, distributed indepen- 
dently of the mx n matrix X and having a gamma density with the Laplace 
transform 

L y (t) = (1 + St)' 13 , 5>0, /3 > (66) 

and let X be distributed as a real rectangular matrix-variate Linnik variable 
as given in Definition 4-1 then the rectangular matrix U = y*X has the 
characteristic function 

<fo(T) = [1 + <5(tr(TS 1 T / S 2 ))f (67) 

As special cases we have the following corollaries: 
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Corollary 4.1 When m = l,Ei = I\ the characteristic function in (67) 
reduces to 

<fo{T) = [1 + <5(T'£ 2 T)?]- /3 (68) 

which is the characteristic function of the multivariable Linnik variable with 
parameters 5,f3, £2 > defined in Lira and Teo (2009). 

Corollary 4.2. When n = 1,^2 = I\ then again (pu(T) is the same as in 
Corollary 4-1 with the parameters 5,(3, Si > 0. 

Corollary 4.3. When m = l,Si = 1,^2 = diag(<7ii, .., a nn ), ajj > 0,j = 
l,...,n then 4>u{T) of (67) is given by 

<j>u(T) = [1 + S(a n h + ... + a nn t n )% ]~f (69) 

Recently Mathai (2010) denned a few matrix random variables in the 
categories of Gaussian Linnik, Gamma Linnik and Kiryakova Linnik. Some 
of these will be given here 

Definition 4.2. The matrix variable U in (67) will be called a real rect- 
angular matrix-variate Gaussian Linnik variable when a = 2 and Gamma 
Linnik variable for the general a. 

By using the multi-index Mittag-Lefner function of Kiryakova (2000) 
one can also define a Gamma-Kiryakova vector or matrix variable. 

Now, let us consider two independently distributed real random vari- 
ables y and X where y has a general Mittag-Lefner density and X has a 
real rectangular matrix-variate Gaussian density as in (62). Let the density 
of y be given by 

..a/3— 1 00 (_-\\k(a\ „.ak 

fM = %r E T(a ; + ak y V>0,S>0,fi>0. (70) 

Then the following are the results given recently in Mathai (2010). 

Theorem 4.2. Let y have the density in (70) and the real rectangular matrix 
X have the density in (62) and let X and y be independently distributed. 
Then U\ = y^X has the characteristic function 



16 



<t>uAT) = [1 + ^MA-'T'B-'T))"}-?. (71) 

Theorem 4.3. Let U2 = y 1 X where y and X are independently distributed 
with y having the density in (70) and X is real rectangular matrix-variate 
Linnik variable with the parameters 7, A, B then U2 has the characteristic 
function 

4>u 2 (T) = [l + ^tr^T'S^T))?]"^, 5 > 0, 7 > 0, > 0, a > 0. (72) 

This U2 will be called a Levy Mittag-Leffler real rectangular matrix-variate 
random variable. 

5. Concluding Remarks 

The Mittag-Leffler functions and their generalizations to vector and ma- 
trix variable cases, discussed in this paper will be useful for investigators 
in disciplines of physical sciences, particularly superstatistics [Beck (2006, 
2008, 2009, 2010)] and nonextensive statistical mechanics [Tsallis (1988, 
2004, 2009)] . The importance of Mittag-Leffler function in physics is steadily 
increasing. It is simply said that deviations of physical phenomena from 
exponential behavior could be governed by physical laws through Mittag- 
Leffler functions (power-law). The paper also continues to show that spe- 
cial functions like Fox H-function and Mittag-Leffler function, represented as 
Mellin-Barnes integrals of the product of gamma functions are suited to rep- 
resent statistics of products and quotients of independent random variables 
[Cottone, Di Paola, and Metzler (2010)]. In this regard, Mathai's (2005) 
pathway model [see also Mathai and Haubold (2010), Mathai, Haubold, 
and Tsallis (2010)] allows to switch from a generalized type-1 beta family 
of functions to a generalized type-2 beta family of functions and further to 
a generalized gamma family of functions. 
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